Symmetry Analysis of the Kohn-Sham Band Structure of Bulk Lithium Fluoride 
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Kohn-Sham orbitals of face-centered cubic lithium fluoride are calculated in prototypical local- 
density approximations to the exchange-correlation functional. The symmetry analysis of these 
Bloch functions in a LCAO basis on a path F-X-W-K-F-L-W through the Brillouin Zone is compiled 
into a list of errata to symmetry labels in the LiF literature, the bulk of which dates back to the 
1970's. 

PACS numbers: 61.50.Ah, 71.20.Ps, 71. 15. Mb 



I. SCOPE 

The crystallographic symmetry of the ground state of 
solid lithium fluoride (LiF) is simple: a face-centered cu- 
bic Bravais lattice with one formula unit in the primitive 
unit cell, known as the rock salt structure, space group 
Fm3m and number 225 in the International Tables 
Yet the symmetry assignments to the lines and points 
of symmetry of (effective) single-particle orbitals (wave- 
functions) populating the reciprocal unit cell vary in the 
literature to a high degree, with no discussion that would 
unify the diverging opinions. This work is to be consid- 
ered an extended erratum to consolidate the graphs of 
band structures of LiF published in the past that are 
incoherent or violate rules of the compatibility tables. 

We start from self-consistent Kohn-Sham (KS) orbitals 
obtained in a Local Density Approximation (LDA) to 
the exchange-correlation (XC) functional of the Density 
Functional Theory (DFT) , and apply the space group op- 
erators to points of low and high symmetry in the Bril- 
louin Zone (BZ), which lead to symmetry labels read from 
the character tables of their little groups. Band structure 
line graphs result as we connect the Kohn-Sham eigen- 
values in accordance with the compatibility tables. 



II. KOHN-SHAM ORBITALS IN THE LDA 
A. GTOFF Bloch Functions 

The basis of this work are the LiF eigenvectors (Bloch 
functions) and eigenvalues (energies) calculated by the 
non-relativistic version of GTOFF Q, H, Q , which solves 
the self-consistent Kohn-Sham equations with an all- 
electron LCAO ansatz for the Bloch functions, and builds 
the atomic orbitals from linear combinations of (Hermite) 
Gaussian Type primitives. We use its "bulk" option 
(translational periodicity of the Kohn-Sham Hamiltonian 
in three directions) and two variants of the LDA, one in 
the next section merely to demonstrate compatibility of 
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the outcome with other, independent calculations, and 
an older one in the main section to make close contact to 
the Zunger-Fremann publication ^ , which, by the num- 
ber of citations received, could be called the "reference." 

Unless noted otherwise, the results are calculated us- 
ing a Kohn-Sham basis with 23 Gaussian Type Orbitals 
GTO's at the Li and 46 GTO's at the F sites (Tabled). 



B. LDA band gap 

The direct band gap Eg at F in the Hedin- 
Lundqvist-type Local Density Approximation (LDA) to 
the exchange-correlation ^ is found to be 

Eg/eY « 8.99 - 3.41A + 1.24A2 (A = a/ao - 7.5939) 

as a function of the lattice parameter 7.37ao < a < 
8.05ao, where a is the edge length of the underlying sim- 
ple cubic unit cell, where a/2 is the shortest distance 
between Li and F, and ag the Bohr radius. (The theoret- 
ical equilibrium lattice constant is smaller than 7.37ao, 
hence underestimates the experimental value, 7.60841ao 
at 25°C 01, by at least 3%. In a recent pseudopotential 
calculation the deviation is 2.7%.) The value is com- 
patible with the 8.7 eV of a pseudo-potential calculation 
with the Ceperley- Alder LDA to the exchange-correlation 
[I at a = 7.62ao, and values of 8.89 eV (8.82 eV) by an 
independent Full-Potential Linearized Augmented Plane- 
Wave program at a = 7.59ao (7.61ao) for the same energy 
functional 

The band gap of the {a — 2/3) calculation is about 
1.0 eV smaller than calculated by Zunger and Freeman 
, which probably indicates that their method of calcu- 
lation starts from parameters of isolated ions and does 
not achieve full self-consistency at the end. (Their pro- 
gram likewise seemed to over-estimate the direct band 
gap of cubic boron nitride [HI by 2 eV compared with 
the consohdated theoretical result of 8.8 

eV, and of diamond .161 by « 0.7 eV compared with other 
theoretical results flTUlsL ITl l2(il| . ) 

A considerably smaller theoretical LiF band gap of 7.65 
eV has been reported by Ching el at. for a = 7.591ao. 
Assessment of this number is not clear; the optical band 
gap of MgO is given as 5.13 eV [23|, whereas an earlier 
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value by what looks like the same methodology is 4.19 
eV H. 



III. SYMMETRY ANALYSIS 
A. Character tabulation 

Point group labels are determined following |23|| . in- 
cluding the Errata by Stokes and Hatch (In addi- 
tion, az ought be added for X in the lattice in p3 . 
Tab. 3.6]. Evarestov and Smirnov interchange Bi and B2 
for Cev in [H Tab. 3.6] relative to [H Tab. 2.2]. This 
is not sir nply a typo but persists if we compare the CW, 
entry in Hi Tab. 4.2] with the Deh table of Flurry |^ 
App. 4]. This does not cause confusion here, as Cqv does 
not enter the subsequent analysis of a cubic space group.) 

The computational procedure is to construct the char- 
acter tables of the ten abstract groups of the lines/points 
of symmetry in the Brillouin Zone according to ,23. p. 
389]. An eleventh entry of the form 

O Gl:{az,0)...:c 

is added to account for points in the 0-plane Tab. 
3.11] that we meet walking from W to K. GTOFF or- 
bitals are rotated to ensure their k-points fall into the 
"irreducible" wedge of the BZ defined in ,23» . Fig. 3.14]. 
For each class of the abstract group of the k-point, a 
character of the GTOFF orbital (Bloch function) (r|i/k) 
of band v is computed as the integral ^^(i^k|{i?|v}|j/k), 
applying one (any) Seitz operator {-R|v} of this class. 
The sum d over up to three degenerate levels deals auto- 
matically with some arbitrariness of the pointing direc- 
tions of the eigenstates/eigenvectors in these cases that 
has been left over from use of a (generic) eigenvalue rou- 
tine in the GTOFF solver. Since GTOFF employs (Her- 
mite) Cartesian Gaussians in the molecular basis, the 
point operators R are adequately represented by Carte- 
sian, matrices: Operation of these on the bases is a multi- 
nomial re-expansion in the a;, y, and z coordinates 
Sec. 1.4]. All operators R are members of the little co- 
group; therefore {i?|v}|j/k) is another Bloch function at 
the same k, and the integrals (matrix elements) become 
standard LCAO overlap integrals. 

This determination of characters associated with a list 
of Seitz operators is numerically stable. Matrix repre- 
sentations of the abstract groups are not needed. It does 
not work blindly in cases of 23, Tab. 5.7] where the 
authors chose generators of the abstract group that are 
not space group operators: consider the example of 
of space group 227: the translation v = has 
been removed from R = C2f, and the overlap integral 
cannot "recover" the missing overall phase factor associ- 
ated with this virtual relocation of the rotated image. (In 
fact, all pair-wise overlap integrals become wrong caused 
by incorrect distances between the Wyckoff positions in 
the "bra" and the "ket.") The alternative is comparison 



of the Bloch function and its rotated image in selected 
regions of the BZ, but this calls for more advanced meth- 
ods of "important sampling" to bypass nodes of the or- 
bitals or dispose regions of small orbital density (in case 
of strongly localized core orbitals). 

The representation of the abstract group that matches 
this list of characters is translated to a point group sym- 
bol by (2^ Tab. 5.8]. Subduction frequencies (compati- 
bility tables) are generated from the point group geneal- 
ogy [2J Fig. 4.1] and [23, (4.3)] for adjacent points on 
the BZ path. 



B. Symmetry assignments and level crossing 

The band structure of Fig. ^ ensues in compliance with 
the compatibility tables ^(5, 27, 28, 29, 30]. Lines are the 
linear (tetrahedral) interpolations/connections obtained 
from the symmetry analysis; so there is no finesse in re- 
gions where perturbation theory proves that there is no 
linear term in the series expansion of the energy as a 
function of k. Point groups labels follow and may be 
converted to the convention of Bouckaert et al. |2^ using 
Tab. ^] The symmetry labels refer to a setting with a 
Li atom at the origin of coordinates, and are not always 
unique. The state near —11 eV at L, for example, con- 
sists of p-orbitals on the F site plus some s-contribution 
on the Li site. This Bloch-type KS orbital is even under 
inversion at this origin. If a F atom is chosen to define 
the origin of coordinates, this same KS orbital is odd 
under inversion at this new origin, and labels have to be 
swapped as follows Tab. 4.7] [31] 32, Fig. 5]: Ai ^Bz, 
A2 ^Bi (at W); A2u ^Mg, Eg ^E„ (at L); A^ B (at 
Q). 

The solid lines in Fig. ^ ^re derived following a no- 
crossing rule for energy dispersions of the same symmetry 
|33| . The rule applies perturbation theory to deduce that 
in almost all cases there is some residual interaction build 
into the Hamiltonian that splits energy levels that close 
up. This argument of [s^, Sec. 1] does not apply to the 
cases of empty orbitals-the simplest demonstration are 
the level crossings in empty lattices (Fig. 4.4.3 in 0|) — 
that are eigenstates of the full/complete (self-consistent) 
1-particle Hamiltonian. Our ground state Hamiltonian 
does not include any mechanism (term) to lift two elec- 
trons to the crossing regions to modify the energetically 
degenerate orbitals. Dashed lines above 25 eV indicate 
some (but not all) alternative choices, if "accidental" de- 
generacies are not ruled out. To decide which of these 
connections are "correct," various approaches are help- 
ful, (i) According to k • p perturbation theory, the over- 
lap (z^kji^'k + q) will be close to 1 in a neighbourhood of 
small q around a k-value, if the correct band v = v' is 
chosen. In practice, however, the density of the mesh in 
the reciprocal lattice hardly is dense enough to have suf- 
ficiently small |q] at hand, (ii) One might compute the 
eigenvalues at an extra intermediate k-point close to the 
expected crossing, (iii) The connectivity of bands does 
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not change when the lattice is set up at a large lattice 
constant (with almost horizontal bands and therefore a 
minimum of accidental crossings), then is "adiabatically" 
contracted to its actual value. With this idea in mind 
we calculate the same crystal at a larger lattice constant 
grown to 110% (Fig.EJ. The prospective accidental cross- 
ing at Q near 28 eV in Fig. ^ is not supported by Fig. 
13 because the state at L is still connected to B2 at 
W near 23 eV in Fig. |2| (This result is confirmed by 
the first approach, the inspection and comparison of the 
dominant basis function of the four states.) In contrast, 
the two accidental crossings near 32 eV at A in Fig. 
seem actually to exist, as a band of Ai symmetry with 
very strong dispersion rapidly moves upwards as a in- 
creases, with its local minimum already above 35 eV in 
Fig. 121 The remaining case nearby 33 eV at A in Fig. 
is not elucidated by Fig. [21 one must increase the lattice 
constant to about 150% to see that the point with E„ 
symmetry at X then falls below the point with Eg sym- 
metry and gets connected with T2(, at F, to indicate that 
there also is accidental crossing of the two lines with E 
symmetry nearby 34 eV on the A line in Fig. ^ 

An additional band of Fj (A2k) symmetry appears in 
calculations with a mixed basis set [s^ and an Aug- 
mented Plane Wave basis (H^, Fig. 1 at 1.9 ryd) which 
is absent in our Figs. ^ and (21 If we augment our LCAO 
orbital basis by a full set of five /-orbitals with exponent 
2.98632 • 10^^ at the F sites, it also emerges from our 
calculations, 18.6 eV above the Fi^c bottom of the con- 
duction band. This basis set completion problem affects 
calculations of dielectric matrices that sum over unoccu- 
pied states [37} . 

The choice of an ATa-type LDA in Fig. ^ and [21 is not 
to promote simplistic density functions, but for immedi- 
ate comparison with the graph in Q and since this could 
serve as a common standard reference of DFT implemen- 
tations. 



C. Discussion 

A tour through the literature in roughly chronological 
order — strictly confined to work that lists eigenvalues or 
band structures — leads to the following comparison: 

• [S^l Fig. 5 is compatible with our results. 

• [i^ The L-labels in Fig. 1 suggest Page and Hygh 
use a setting with Li in the center of coordinates: 
Wi near -0.45 ryd should read Wj. 

• [s^ The results in Tab. 1 and Fig. 3 are compati- 
ble with ours. Chaney et al. did not include d-type 
basis functions, which means the T2£, point near 15 
eV in our Fig. ^ and bands that are attached to 
it disappear. In consequence they neither show a 
Fjsj, nor a in Tab. 1. Our Fig. [3 illustrates this 
case. The d-orbitals are obviously a key ingredient 



to construct an even function at F in the conduc- 
tion band, with otherwise only s-hybrids remaining 
which have to stay orthogonal to the three s-type 
core orbitals. 

• [S^ Figs. 1, 3 and 6: the upper should read 
Li. (The compatibility relations would allow the 
connection F1-A1-L2.) The equivalent correction 
in Tab. VIII is to exchange L3 — > with L3 Li 
in the first column. 

• ^1 The results in Tab. 1 are compatible with 
those here in the same sense as discussed for 38]. 

• Figs. 1 and 2: the upper should read Li. 
Tab. 1: replace by Lie. 

• Figs. 3 and 4 are compatible with our analysis. 

• Q Fig- 1: the 1 at X should read 4'. A comment 
on the 3 at X and L is difficult because the order of 
their Hartree-Fock energy levels might be switched 
in comparison to our DFT results. 

• m The label X3 in Tab. IV should read X^ where 
they refer to their own results, as their basis does 
not include d-orbitals, see the discussion of [ssj'). 

• HBl Figs. 2 and 3: the uppermost X4 should read 
Xi. (The compatibility relations would allow the 
connection F15-A1-X4.) The uppermost L3 should 
read L3. The middle Lj near 10 eV should read Li. 

• 113 In the LiF part of Fig. 2, the 1 at X should 
read 4'. (The compatibility relations would allow 
the connection Fi-Ai-Xi.) 

• ^3 This tight-binding analysis of the valence bands 
confirms our results. 

• ^3 The lower L3 in Fig. 8 should be replaced by 
L'2, and the Xi by X^. 

• The labels at L should read (bottom- up) Li, Lj, 
L3, Li, L2, L3, L2, and L3 in Fig. 2. All three cases 
of L3 and L3 should split towards Q. There should 
be no crossing at Q near 7 eV. The A-lines near 
12 eV should twist to connect X3 to F25 by A2 and 
Xi to Fi5 by Al. The labels ^2 and S3 near 17 
eV should probably be exchanged. In Tab. 1, all 
X4, X5 and L2, also L3„ should be primed. Xi^c 
in Tab. 1 should read Xs^c where it refers to their 
own calculations, which follows from measuring the 
vertical distance in Fig. 2 and comparison with the 
value of 0.710 hartree given in the table. (See again 
the comment on [13 •) 

• [sj Fig. 3 is compatible with our labeling scheme. 

• [H^l In Fig. 2, the upper Lj should read Li. 

• lU Fig. 1 is reproduced from hence the same 
comments hold. 
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• [3 The following LDA lines in Fig. 1 should cross: 
Z lines near 41 eV; O lines near 41 eV; S lines near 
30 eV and 39 eV; S lines close to K near 35 eV. 



tenance of the GTOFF code was an indispensable foun- 
dation of the present analysis. 
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TABLE L Exponentials and contraction coefficients of GTO's 
centered at F and Li, derived from Partridge's 13s9p and 
15sl0p sets lillHl. 
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TABLE IL The interpretation of the point group symbols 

of ;2l 

in terms of those used in l29l . obtained by searching 
for matching lines in the character tables. There is only one 
table for mm2 in I2M p. 58] and only one common table for 
Z and K in |29l . Tab. VI], but the ax and a^, operations are in 
different classes in l2g| . Tab. VI], which leads to two separate 
columns for Z and K here. 




FIG. 1: KS band-structure of LiF at o = 7.5939ao = 4.01852 
A with the Xa form of the XC energy (a = 2/3). Crosses 
mark energies in the BZ that were calculated by GTOFF on 
a grid of 16 x 16 x 16 points in the reciprocal unit cell. The 
direct gap at F is 8.79 eV. 



FIG. 2: KS band-structure of LiF as in Fig.[T]for a 10% larger 
lattice constant. As the lattice constant increases, the bottom 
of the conduction band at F becomes a local maximum with 
satelhte minima at the adjacent A, E and A lines. 
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FIG. 3: KS band-structure of LiF as in Fig. but the basis 
of Tab. H] replaced by a smaller basis of Li lls4p (contracted 
41111111/1111) and F 13s8p (4111111111/2111111) to illus- 
trate the role of F d-orbitals. 



